NASA TN D-1447 



Nb3~/036f 

COBS I] 

NASA TN D-1447 



TECHNICAL NOTE 


D-1447 


THE MAGNETIC FIELD OF A MODEL RADIATION BELT 
NUMERICALLY COMPUTED 


Syun-Ichi Akasofu 
University of Alaska 
College, Alaska 

Joseph C. Cain 
Goddard Space Flight Center 
Greenbelt, Maryland 

Sydney Chapman 
University of Alaska 
College, Alaska 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON November 1962 



THE MAGNETIC FIELD OF A MODEL RADIATION BELT 
NUMERICALLY COMPUTED 


by 

Syun-Ichi Akasofu 

University of Alaska 

Joseph C. Cain 
Goddard Space Flight Center 

Sydney Chapman 

University of Alaska 


SUMMARY 

The magnetic field of a ring-current belt symmetrically encircling 
the earth is calculated numerically, to a first approximation, for a par- 
ticular model belt discussed elsewhere by Akasofu and Chapman. They 
calculated the first approximation to the field only for points in the 
equatorial plane. The whole distribution of the field is discussed here, 
and for a particular intensity of the belt the calculation is carried to a 
second approximation. The field of the ring current is shown to be nearly 
uniform over the earth 1 s surface, although the diamagnetism of the belt 
produces an insignificant irregularity in the field disturbance in auroral 
latitudes. In the model belt considered here is one whose center line is 
connected with the auroral zone by the dipole field lines of force. Although 
its particles have not yet been observed directly, its existence is suggested 
by the quiet day anomaly of the satellite-observed magnetic field at six 
earth radii. The known Van Allen radiation belts probably contribute 
little to the ring current field, but during magnetic storms the radius of 
the belt mainly responsible for this field is probably less than this dis- 
tance. 
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THE MAGNETIC FIELD OF A MODEL RADIATION BELT 
NUMERICALLY COMPUTED* 

by 

Syun-Ichi Akasofu 
University of Alaska 

Joseph C. Cain 
Goddard Space Flight Center 

Sydney Chapman* 

University of Alaska 

INTRODUCTION 

Space observations by Van Allen and Frank (Reference 1) and by Vernov, Chudakov, 
et al. (Reference 2) have shown that the earth is encircled by radiation belts which are 
composed of charged particles trapped in the magnetic field f in the space around the earth 
(Reference 3). 

A group of particles trapped in the earth’s magnetic field produces a magnetic field AF 
by their motions in this field. The total field vector B is the vector sum of the intensities 
B e and AB of the earth’s magnetic field f e and the field AF of the belt; thus 

H = H E + AH , (!) 

which may be expressed symbolically by 

F = F e + AF . ( 2 ) 

The total field F and the precise distribution of the particles in the belt are intercon- 
nected in a somewhat complex way. Successive approximation provides one method of 
investigating this relationship. To apply this method to a model belt whose particle content 
and general distribution are known or ’’assumed”, the following steps are taken: 

1. Calculate the particle distribution and the magnetic field of the belt as if the total 
field present were simply the earth’s field F E . Let the belt field thus calculated be 
denoted by AjF. Let the field F E + A 2 f be denoted by F t . 

2. Calculate the particle distribution and the magnetic field of the belt as if the total 
field present were F r Let the belt field thus calculated be denoted by A 2 f, and let 
F 2 denote the field f e + A 2 f. 

•This report has been published in substantially the same form in J. Geopbys. Res. 66(12):4013-402 6, December 1961- 
tAlso engaged at the High Altitude Observatory, University of Colorado, in a program of research sponsored by the National 
Bureau of Standards and the Air Force Geophysical Research Directorate. 


1 


2 


By continuing this process, a series of values of the belt field A n F and of the combined 
field F n are obtained: 


If F n approaches a limiting value, this is expected to be the true field F given by Equation 1. 

Akasofu and Chapman (Reference 4) have discussed the field A 2 f analytically for sev- 
eral model belts. They treated f e as the field of a point dipole of vector moment M 
(8.31 x 10 25 gauss cm 3 ) situated at a point 0 inside the earth (centered or eccentric). In 
this case the belt and its field AF, hence the total field F, have symmetry about the dipole 
axis and equatorial plane. Also they made numerical calculations of several properties 
for a particular assumed model belt called v 3 , including the intensity A,h of the field A 2 f, 
for points in the equatorial plane of the dipole M. The results obtained were compared with 
the Explorer VI (1959 6) magnetic observations (Reference 5). 

Their model belt involves a parameter n 0 E which, for the same distribution of the 
belt, determined its density everywhere. The field A 2 f involves n 0 E as a factor; that is, 
its intensity everywhere is proportional to n 0 E . For the particular value n 0 E =150 the 
calculated field F, seemed to resemble the field observed by the satellite Explorer VI in 
the region up to about 7a (where a = 6.37 x 10® cm is the earth’s radius). However, in 
the region between 5a and 7a, the observed field showed a notable proportionate departure 
from the earth's field f e . Hence the calculated field A 2 f may not be a good approximation 
to the actual field AF of the belt. Akasofu and Chapman (Reference 4) stated that "further 
work on this problem is necessary to obtain a self-consistent system of current and field." 

In this paper their numerical calculations for the model belt v 3 are extended to show 
the general character of the field A 2 f in the whole space around the earth. These results 
are compared with the changes observed at the earth's surface during magnetic storms. 
Also, the second approximation A 2 f to the field of this model belt has been calculated, for 
n 0 E = 90, in a region that includes the densest part of the belt. The difference between 
AjF and A 2 f suggests that, for this value of n 0 E, the field A 2 f is a good approximation toAF. 

This problem is being explored further, particularly by calculating the potential of AF 
(or at least of A 2 f) in terms of zonal harmonic functions. Alternative values of the belt ra- 
dius and thickness, and also of a parameter of the pitch angle distribution are being 
considered. 

GEOMETRICAL CONSIDERATIONS 

The point O at which the dipole M is located is taken as the origin of coordinate systems 
r, 6 , X and x, z , X , relative to a polar axis Oz antiparallel to M; z is measured parallel and 
x perpendicular to Oz; X increases eastward, z northward. In a magnetic field which, like 
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F e andF, is symmetrical with respect to Oz and the equatorial plane z = 0, it is convenient 
to substitute, for the coordinates of any pointP, the parameters r e , s w where r e de- 
notes the distance from O to the point P e (where the line of force through P crosses the 
plane z = 0) and s x denotes the length of the line of force from P e toP (positive northward). 
Let h, j , k denote a right-handed orthogonal triad of unit vectors, with origin at P, •> being 
the direction of H at P, and k along the eastward normal to the meridian plane through P; 
j = kxh . An element of length d« drawn fromP is 

ds = h dSj + jh 2 dr e + kh 3 dX. . (4) 


Areal elements at P normal to k, j, k may be denoted by dS, , dS 2 , dS 3 : 


ds l 

= h 2 h 3 dr e d\ . 


dS 2 

- h 3 ds 1 dX , 

( 5 ) 

dS 3 

= h 2 dr e dSj . 



Let H e , H denote the (scalar) magnetic intensities at P e and P. The constancy of the strength 
HdSj of a tube of magnetic force along its length gives the relation 


h - ^ 
h 2 “ h,H 


h 3 - r sin 6 . 


( 6 ) 


Akasofu and Chapman (Reference 4) gave the special forms of h 2 , r e , H, H e , and the radius 
of curvature R c of the line of force at P, for the field F E of the dipole M, in terms of r 
and 6 . These factors are used in evaluating AjF. In their Equation 4 hj (=“h x d^) was 
used instead of ds 1? where 4 > denoted the geomagnetic latitude. 

In addition to the general relations of Equations 4 to 6 we also note that 

v = h A, + j + k "h faT • (7) 


and 


BX(HV)H 



( 8 ) 


Equation 8 refers to a field like F, symmetrical about Oz, whose lines of force lie in me- 
ridian half -planes. 
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THE MODEL RADIATION BELT 

The belt can be specified at each point P (with position vector r relative to o) by the 
function n(r, w) expressing the distribution of particle velocities w . There is such a 
function for each kind of particle in the belt; n(r, w) dw denotes the number density of 
these particles, at r, whose velocities lie in the small range dw about the value w. By 
symmetry, n is a function of s 1? r e , and w; hence it can be written n(s 1 , r e , w) . 

The angle between w and H at P is called the pitch angle 4 > . The distribution of the 
azimuth of w around H is isotropic except possibly at singular points of the field F . Hence, 
for the particles of this velocity group, the velocity distribution depends only on It may 
be specified by a function F(<£) such that F(<£) d^ denotes the fraction of all the particles of 
this group whose pitch angles lie between <f> and c £ + d<£. Thus 

| F(0)d0 = 1 ; (9) 

Jo 

and the number density at P of the particles (of any particular kind) with speeds between 
w and w + dw can be denoted by n(Sj, r e , w)dw. Clearly 

j n(r , w ) dw = 2 rr j* jn( s t , r e , w) F( 0 ) w 2 dw d<£ 

" 277 f n(s 1( r e , w) w 2 dw . (10) 


In general, as the particles move in the field F, their pitch angle distribution and 
number density vary from point to point. Parker (Reference 6) showed that the pitch angle 
distribution given by 


F(<£) - A(a.) sin a + 1 cfi , 


where 


(ID 


A(a) 


2 a + 


r<a + 2) 



which satisfies Equation 9, has the following remarkable properties: 

1. If F ( 0 ) has this form at one point along a tube of force, this is its form for the 
same value of a all along this tube; and 
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2. The number density varies as H' a/2 along the tube, and we may write 

/ H eV /2 

n ( Sl> r e ,w) = n(r e , (12) 


where n(r e , w) dw denotes the number density at P ( of the particles with speeds 
between w and w +dw. Clearly the functional symbol n in Equation 12 has different 
meanings in the two cases. 

Following Akasofu and Chapman (Reference 4) we consider a model belt in which F(d>) 
has the form of Equation 11, with the same value o throughout the whole field. Further, 
the density distribution in the equatorial plane z Q 0 is 

n(r e , w) = n Q (w) exp [-q 2 (r e ~ r eQ ) 2 ] . (13) 


Thus in this plane the radial distribution of the density is Gaussian for each particle speed. 
The factor q determines the radial spread of the belt. Conceivably q might be a function 
of w, but in this paper we consider only particles of a particular speed w or speed range w to 
w + dw ; thus the results should be integrated with respect to w if the distribution functions 
of the speed w(or energy E = 1/2 mw 2 ) are known for the particles present. 

The model belt of Reference 4 is specified here by Equations 11 and 13 with 



These parameters are chosen merely to illustrate the numerical computation, and the same 
procedure can be applied to any set of parameters. In this case, the dipole line of force 
that crosses the plane z = o at 6a meets the earth 1 s surface in geomagnetic colatitude 
0 = 24° . The choice 6a for r n , the radius of the center line of the belt distribution, is 
made because aurorae often appear during moderate magnetic storms at about that colati- 
tude. The value a = - 1/2 corresponds to a pitch angle distribution that has more small 
pitch angles than would correspond to an isotropic distribution of the direction of w. The 
chosen value of q corresponds to a reduction of density to a tenth of the maximum value 
(also denoted by n 0 (w)dw or simply by n 0 ), at a radial distance a on either side of the cen- 
ter line of the belt, that is, at 5a and 7a. 
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THE CURRENT DISTRIBUTION CORRESPONDING TO AJ 

Let i denote the contribution to the equivalent electric current intensity in the belt 
due to particles of any one kind with speeds between w and w +dw. Let i t denote the cur- 
rent intensity for the first approximation. Following Parker (Reference 6), Akasofu and 
Chapman expressed i x for a belt with the pitch angle distribution of Equation 11 substan- 
tially as follows (Reference 4, Equation 51): 

mcw 2 dw/ H e \° /2 [ , [l - 6B(a)~| dn(r e , W ) 2^) 1 

" H \hJ \ n ( r e’ w ) R c dr e h 2 J 


where 


(a+2) 

B ( a ) “ 4(a + 3) 


(16) 


Equation 15 is valid for the general field F, where the current flows symmetrically round 
the axis Oz ; such a current distribution is called a ring current. 

Akasofu and Chapman calculated the distribution of i for the model belt specified in 
Reference 4, section 3* when the total field present is taken to be the dipole field F E . 

In this case their i will be denoted by i 1 ; then at any point P 0 on the center line of the 
belt (z = 0, r e ~ r e0 )} the value of i 1 is denoted by 


i j 0 - 1 . 56 x 10 3 mw 2 n Q (w) dw esu 

= 1.66 x 10“ 15 n Q (w) E dw amperes 



J 


(17) 


where E is the kinetic energy of the particle expressed in kev, 

E = 3.12 x 10 8 mw 2 . (18) 

Akasofu and Chapman (Reference 4, Figure 3c) illustrated the current distribution by 
a set of current isolines, for the current intensities i x /i 1 0 = 12, 10, 8, 6, 4, 2, -2, “4, -6 . 
Some of these lines (10, ±6, ±2) are reproduced here in Figure 1; the positive values corre- 
spond to westward current, negative to eastward, and the eastward isolines are dashed lines 
in the Figure. (Their Figure 3c shows by dot-and-dash lines the dipole lines of force 
through the points r e = 5a, 6a, 7a .) The current intensity changes sign near P 0 ; slightly 

*Note that a factor (1 + sin 2 4>) should be added on the right of their Equation 53, that a factor cos 4> is omitted from Equation 
29, and that in Equation 51 the index of the factor H e /H should be a/2. The numerical factor 1.66 x 10 1 in Equation 17 
was written as 1.61 x 10 “ 1S . These are errors only in the reproduction of the original text; the calculations were made 
from the correct formulas. 
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nearer to 0 ; the maximum westward and eastward current intensities slightly exceed 
12 i 10 and 6 i 10 . Also they illustrated separately (Figures 3a, 3b) the distribution of parts 
of i given in Equation 15; the distribution of the parts of i that are contributed by the 
diamagnetism and drift currents of the belt were to be shown later. The main contribution 
to the field comes from the diamagnetism. 



Figure 1— The field vectors of A^, the first approximation to the magnetic 
field of the model belt, in a meridian plane; and isol ines of the equivalent 
current intensity in the belt (solid lines indicate westward, dashed lines 
eastward current). The vector scale of force, and the unit in which the 
current intensity is expressed, are proportional to the energy density n Q E 
(kev/cm 3 ) at the center line of the belt, at 6 earth radii from the earth's 
center. 


THE MAGNETIC FIELD A,F 

Akasofu and Chapman gave expressions (Reference 4, Equations 76 and 77) for the x 
and z components of the field A X F , and calculated A^ numerically for points in the equa- 
torial plane z = o , where A j h x = o (see their Figure 6). Using the same formulas, and 
calculated values of i x at a network of points P in a meridian half-plane, the present authors 
evaluated A 1 H x and A l H z at a large number of points P' in this half-plane. Figure 1 shows 
some of the results by the vectors A x h for the points p' disposed along radii in latitudes 
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4> 1 from 0° to 70° at 10° intervals, mainly at intervals (a) of radial distance. Figure 2 
illustrates these and other results in a different way. The lower portions of the figures are 
graphs of the x and z components of the A,F field on a linear scale along radii in latitude 
4> 1 from 0° to 50°, up to r = 8a , for a model belt such that n 0 E = 150, where 

n o = n o( w > dw • (19) 

This might correspond, for example, to a maximum number density 1.5/cm 3 of 100-kev 
particles, or 1/cm 3 of 150-kev particles, at points P 0 on the center line of the belt. The 
upper parts of Figure 2 are graphs, on a logarithmic scale, of the intensity of the dipole 
fields F e and F, . Figures 2b and 2c (upper) also show graphs of the x and z components 
of the field F 1 . 

For this value of n 0 E it is clear from Figure 2 that in certain regions, mainly near 
the center line of the belt, AjF/f e is not small. Hence for this belt intensity it is inadequate 
to ignore the distortion of the field F E by A,F . However the graphs in the lower parts of 
these figures are valid for smaller values of n 0 E, with scales of force reduced proportion- 
ately. 

Figures 1 and 2a also indicate that within a radius of about 2a the field A t F is nearly 
uniform. Beyond this \F becomes notably nonuniform, and shows considerable curl in 
the region of most intense current, centered about 6.5a in the equatorial plane. The field 
differs considerably from that of a toroidal current. Outside the belt the field A t F tends, 
with increasing distance, to similarity with the dipole field, as general theory would indi- 
cate. 

Figure 2 may be useful for comparing the theory of the field of the model belt here 
considered with observations of the actual magnetic field around the earth. Clearly the 
greatest magnetic influence of the belt should be explored by satellites moving nearly in 
the equatorial plane. For comparison with satellite magnetic observations, tables of AjF 
have been calculated for many more points than are illustrated in the figures of this paper. 

Figures 3 and 4 show some lines of force of the combined field F, = f e + A,f. These 
lines were calculated approximately, taking steps of length o. la , for the values n 0 E = 90, 
150. These lines correspond, for example, to distributions of 150-kev particles, of maxi- 
mum number densities 1/ cm 3 and 0.6/cm 3 . Figures 3 and 4 show typical lines of force of 
the fields F t (solid lines) and f e (dashed lines). The distortion of the F t lines is much 
greater for n 0 E =150 than for n 0 E = 90. Within the center line of the belt, the ring 
current reduces the magnetic flux across the equatorial plane. Outside this line the flux 
is increased by an amount depending on the magnetic moment of the belt. Akasofu and 
Chapman gave a general formula for the magnetic moment of a belt of the type here con- 
sidered (see Reference 4, section 4.4). 


(•K>d|p) j 



the disturbance vectors in the x and z directions. 






DISTANCE r. (earth radii) DISTANCE r. (earth radii) 











n 0 E = 90 


Figure 3—The lines of force of F , the first approximation to the combined field of the geomag 
netic dipole and the belt. The dashed lines are lines of force of the dipole field. 


n o E=150 



: igure 4 — The lines of force of F , the first approximation to the combined field of the geomag 
netic dipole and the belt. The dashed lines are lines of force of the dipole field. 
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THE RING-CURRENT MAGNETIC DISTURBANCE OF 
THE EARTH'S SURFACE FIELD 

The ring -cur rent field at the earth’s surface is nearly uniform and in the direction 
of M , the dipole moment, as is shown particularly by Figures 1 and 2a. Thus it decreases 
the horizontal component of the field and increases the numerical intensity of the vertical 
component. Changes in the field are ob- 
served during the main phase of magnetic 
storms. As the ring- current grows during 
the early hours of the storm, its changing 
field induces electric currents within the 
earth. They increase the change in the hori- 
zontal component and reduce that of the 
vertical component. 

In Figure 5 the ring-current change of 
the horizontal component (denoted by u) at 
the earth’s surface, is given as a function 
of latitude up to 70°, for n 0 E = 150. On 
the scale shown, the curve is indistinguish- 
able from a sine curve corresponding to a 
uniform field suggested by Figure 1. The 
diagram also shows the graph of the verti- 
cal component v of the ring-current field 
at the earth’s surface for n 0 E = 150. The 
curve is almost indistinguishable from the 
cosine curve corresponding to the sine 
curve for h, except near 65° latitude. 

However, calculations of the field components of AF were made for many points within 
and near the ring-current distribution (Figure 1), and they revealed a small irregularity of 
the field f in the current region, extending even to the earth’s surface. This is to be attri- 
buted mainly to the diamagnetism of the belt. Figure 6 illustrates this feature in closer 
detail. It shows the variations of the total intensity A^, and of its horizontal and vertical 
components H and V at the earth’s surface, between latitudes 58° and 70°. The irregularity 
affects V most, but its magnitude for n 0 E = 150 does not exceed 5 gammas (Reference 7, 
p. 2250). This model belt would reduce the earth’s horizontal magnetic intensity at the 
equator by about 50 gammas; the induced earth currents would increase this to about 90 
gammas (and would decrease the change in the vertical component). The irregularity of 
the vertical component in the auroral zone will not exceed 35 gammas. The observed de- 
crease of the horizontal component at the maximum phase of the greatest magnetic storms 



Figure 5— The variations of the horizontal and 
vertical components (H, V) of the field of the 
belt, at the earth's surface, from the equator to 
latitude 70', to the first approximation for 
n 0 E - 150 
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is of the order of 600 gammas. The actual 
changes in the vertical component are often 
as great as 1000 gammas, and sometimes 
much more. In those latitudes the field of 
the belt, mainly diamagnetic, is insignifi- 
cant by comparison, as was recognized by 
Dessler and Parker (Reference 7). This 
conclusion contrasts with the diamagnetic 
explanation of polar magnetic disturbance 
suggested long ago by Hulburt (Reference 
8) and criticized by Chapman (Reference 9, 
section 10). 

THE SECOND APPROXIMATION (A*F) 

TO THE FIELD OF THE RING CURRENT 

For the values of n 0 E considered in the 
two preceding sections the field of the ring 
current materially modified the dipole field 
in the regions of greatest current intensity. 
Hence our calculation of AjF, the first ap- 
proximation to the ring-current field, needs 
correction for comparison with observed 
magnetic changes during storms. 

The current intensity has been recalcu- 
lated at a network of points P m n (Figure 7), 
for the model belt specified by Equations 
11, 13, and 14, in the presence of the field 
Fj = F e + AjF . This calculation has been 
made for the case n 0 E = 90 . The point 
P lies on the line of force for which 

m , n 

r e = 4.8a + m x 0. la and is at the distance 
s x along this line from the plane z = 0; s 2 
is given by = 0. la + n * 0. 2a. The values 
of m ranged from 0 to 27 (i.e., 28 lines of 
force), and the values of n ranged from 0 
to 27 for m = 0 and from 0 to 45 for m = 27. 

The calculated current intensity at P m n 
is denoted by (i 2 ) m n , where the subscript 



Figure 6— Enlarged version, for the latitude range 
58° to 70° (near the tip of the model belt), of the 
H and V curves of Figure 5; and th e variatio n of 
A t F, the corresponding total vector VTP _ + _ V T (the 
dotted line shows the interpolated curve for AjF, 
to indicate the si ight irregularity due to the belt). 
Note the three different scales for H, V, and AjF. 



i i 

I i 

\ i 


0. la 


Figure 7— The subdivision of the meridian plane 
used in computing A 2 F, the second approximation 
to the field of the belt; the curved lines are lines 
of force of the field Fj (see Figures 3 and 4). 
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2 indicates that the value is a second approximation to i , corresponding to motion of the 
particles in the field Fj . The values of i 2 were calculated from Equation 15, using the 
values of H, h 2 , R c numerically computed for the field F 1 and its lines of force. 


The values of i 2 at the points P m n having been obtained, the magnetic field A 2 f was 
computed at another network of points n , by using formulas 76 and 77 of Reference 4. 
The current distribution was divided into filaments (Figure 7) centered at the point P 
and of cross sections h 2 dr e ds x ; the value (i 2 ) m n was taken to represent the mean cur- 
rent over this section. The number of points P m n in the two (northern and southern) halves 
of the meridian section of the belt exceeded 2000. The field A 2 f was computed for points 
P' m n along radii in latitudes 4> x - 0°, 10°, 20°, at radial distances a to 8a. The radial 
interval was 0.2a from r e = 4.4a to r e = 7.6a. These results are embodied in numerical 
tables (not reproduced here). 


Figure 8 illustrates the variation 
of A 2 f and F 2 in the equatorial plane 
in the same manner as Figure 2a does 
for AjF and F x ; Figure 2 refers to 
n 0 E = 150, and Figure 8 to n 0 E = 90. 

The second approximation corrects 
\F in this plane by a maximum 
amount of 15 gammas, at a point near 
P 0 . The change from AjF to A 2 F | 

scarcely alters the location (just | 

within the center line of the belt) 
where A t F is most intense, but it \ 

deepens this minimum by about 15 
gammas, from 75 to 90. Except in 
this region, which produces a mini- 
mum of Fj and F 2 near P 0 , the 
second approximation makes a cor- 
rection of less than 10 gammas; and 
at the earth’s surface the correction 
does not exceed 3 gammas. It seems 
likely that further approximations 
would produce little additional cor- 
rection, so that A 2 f and F 2 nearly 
equal AF and F . 

The dip in the fields F x and F 2 
between 5a and 7a is shown on a 
linear scale by an inset in Figure 8. 



Figure 8-Similar to Figure 2, but for n Q E = 90, and with 
the addition of the second approximation toAF (below) 
and to the combined field F (above). The broken line 
gives the first approximation, the solid line the second. 
Top right: the variation of F x (broken line) and F 2 (solid 
line) on a linear scale, between 5a and 7a. 
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Near the center line of the belt, for n 0 E - 90, the estimated intensity of the true field 
F in the equatorial plane will be reduced to about 55 gammas. 

For n 0 E = 150, the first approximation to the combined field is shown in Figure 2a; it 
reaches a minimum value of about 24 gammas. Because in the second approximation, the 
minimum value is smaller than that for the first approximation, it may be expected that 
the above value of 24 gammas also becomes smaller. 

If n 0 E is great enough to reduce the field F to a low value near p Q , the guiding center 
approximation (Reference 10) used in obtaining Equation 15, may cease to be valid near 
P 0 . The paths of the particles there would be extremely complicated and the current inten- 
sity calculation must be modified. However, the dip near p Q is mainly produced by the 
strong eastward current and the westward current, centered respectively at about 5.5a 
and 6.5a, where AF is relatively small and where the guiding center approximation will re- 
main approximately valid. 

In order to demonstrate this, A 2 f was recalculated for n 0 E = 90, neglecting the con- 
tribution from the current in a region around P 0 , limited by the lines of force r e = 5.9a 
and r e = 6.1a, and by the radii between latitudes ±10°. At 6a, the value of A 2 f thus calcu- 
lated was found to be 83.9 gammas; this may be compared with A 2 f = 90.1 gammas for the 
complete belt. Thus, more than 90 percent of the field decrease near p Q is produced by 
the particles outside the above limits. 

There is some doubt about what will happen if the value of n 0 E is further increased. 
Beyond a certain value of n 0 E the field may be reversed near p Q . It is interesting to note 
that such reversals have been observed in the thermonuclear reaction vessels (mirror 
machines) in which a very hot plasma is confined by a strong magnetic field (Reference 11). 
Akasofu and Chapman (Reference 12) have attributed much importance to such a reversal 
in the formation of auroras. 

THE RING-CURRENT FIELD BEYOND THE BELT 

Akasofu and Chapman calculated the first approximation a 2 m to the magnetic moment 
of the model belt considered here (Reference 4, section 4.4). Their value was 1.17 x 10 23 
n 0 E gauss cm 3 . Hence the equatorial field intensity F x beyond the belt will tend asymptoti- 
cally with increasing distance to the value (M + A x M)/r e 3 . Table 1 gives the values of 
M/r e 3 , (M + A 1 M)/r e 3 , and the calculated values of F x , from 8a to 15a at intervals of a for 
the case n 0 E = 150. The values are expressed in gammas. Table 1 shows that at 15a the 
field approximates closely that of a dipole m + A x m . 

For n 0 E = 90, equatorial values of F 2 have been calculated up to r c = 8a. The values 
of Fj and f 2 at 8a in this case are 80.0 and 81.7; at this distance (M x + A x M)/r e 3 for n 0 E = 90 
is 70.4. 


Table 1 

The ring- current field beyond the belt 
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Magnetic 

Parameter 

Calculated Value* (gammas) 

r e = 8a 

r e = 9a 

r, = 10a 

r e = 11a 

r e = 12a 

r e = 13a 

r e = 14a 

r e ~ 15a 

M/r e 3 

62.5 

43.9 

32.0 

24.0 

18.5 

14.6 

11.7 

9.5 

(M + AjM)/r e 3 

75.7 

53.2 

38.8 

29.1 

22.4 

17.6 

14.1 

11.5 

F i 

91.6 

61.1 

i i 

43.2 

31.8 

i 

24.2 

18.8 

15.0 

12.1 


•Reference 4 


THE NONSTEADY STATE 

During a magnetic storm there are great changes in the number and distribution of 
the energetic particles in the region beyond the inner radiation belt. Akasofu and Chapman 
(Reference 4) have discussed these changes, as observed from satellites (Reference 13) ? 
in connection with the magnetic variations during particular storms. 

The formulas used to determine the field of the ring current of a steady radiation belt 
do not suffice for a changing belt. The corresponding changes in the field f introduce an 
electric field e such that 

dH 1 

Tt ~ ~~ curl E • 

Assuming axial symmetry, and using cylindrical coordinates, 

3l - (±\ 8(rE ^ 

dt ~ \ cr / 

for the electric intensity E^ in the direction k . This electric field adds a term v equal to 
(EXH) /cH 2 or E^j cH z to the velocity of the particles, and a term ( nmc H 2 ) H Xdv / dt to the 
current intensity i. The velocity v is directed outward or inward, in the meridian plane; 
it tends to distort the density distribution. A complete discussion of the ring current must 
take account of this electric field and of various loss processes in the belt. 


ADDENDUM 

Further computations of the magnetic field have been made for measured particle dis- 
tributions during quiet (Reference 14) and disturbed (Reference 15) magnetic conditions. 
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Computations have also been made (Reference 16) for different types of belts. Further 
discussion of the methods of calculation of the magnetic field of trapped particles is given 
in References 17 and 18. 
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